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ABSTRACT 

We present a detailed study of the real and integrals-of-motion space distributions 
of a disrupting satellite obtained from a fully self-consistent high-resolution cosmolog- 
ical simulation of a galaxy cluster. The satellite has been rc-simulatcd using various 
analytical halo potentials and we find that its debris appears as a coherent structure 
in integrals-of-motion space in all models ("live" and analytical potential) although 
the distribution is significantly smeared for the live host halo. The primary mecha- 
nism for the dispersion is the mass growth of the host. However, when quantitatively 
comparing the effects of "live" and time- varying host potentials we conclude that not 
all of the dispersion can be accounted for by the steady growth of the host's mass. 
We ascribe the remaining discrepancies to additional effects in the "live" halo such as 
non-sphericity of the host and interactions with other satellites which have not been 
modeled analytically. 

Key words: galaxies: clusters - galaxies: formation - galaxies: evolution - n-body 
simulations 



1 INTRODUCTION 

With the d iscovery of st e llar s tr eams within the Milky 
Way (e.g [Helmi et al. I Il999t IChiba fc Beers! I 200C : 



Ibata et al. 1(2001 iBrook et al. 1 120031: lYannv et al. I 



2003 



Navarro. Helmi fc Freeman I 120041: IMaiewski et al. I 2004 ) 
and M31 (e.g. Ilbata et al. I l200ll iMcConnachie et al. I 
20041) and streams and shells in clusters 



Trentham & Mobasher 
Calcaneo-Roldan et al. I 



Il99gl : 

2000; 



Gregg & West 



Fcldmcicr ct al. 




[ Helmi fc de Zeeuw Il200ot lHarding et al. Il200lt Ilbata et al. I 
120021 iBekki et al. Il200a IMaiewski et al. 112004) . Such sim- 
ulations have prov ided great insight into the forma- 
tion of stellar (e.g. Ilbat a et al. I l2002t) and gase ous (e.g. 
lYoshizawa fc Noguchi l l2003l:IConnors* et al. 1 l2004l) streams, 
and helped to constrain the shape of the Mi lky Way's halo 
(e.g. Ijohnston et al. 1999; Ibata et al. 
ture content (e.g. I Johnston et al. 1 12002') 



2001) and substruc- 



they have become a standard fixture in our understanding 
of galaxy and cluster formation. These streams provide 
important observational support for the hierarchical build- 
up of galaxies and clusters and the A-dominated cold 
dark matter (ACDM) paradigm. Future observational 
experiments such as RAVE 1 and GAIA 2 are designed to 
shed further light on the importance of streams in galaxy 
formation. Thus to better interpret this observational data 
coming online in the near future it is necessary to have 
clear theoretical understanding of the stellar streams left 
behind by dissolving satellites. 

The most common approach to this has been to simu- 
late the disruption of individual satellite galaxies in static 
analytical potentials representa tive of the distribution of 
dark matter (DM) halos (e.g. iKoiima fc Noguchi I Il997t 



http : //astronomy . swin . edu . au /RAVE/ 



http : //astro . e stec . esa . nl/GAIA/ | 



The ACDM structure formation scenario predicts that 
small objects form first and subsequently merge to form en- 
tities thus DM halos are never "at rest" but always be in 
the process of accreting material from its vicinity and hence 
grow in mass. This also extends to the Milky Way, although 
it e ncountered its last major merg er some 10 Gyrs ago 
(e.g. iGilmore. Wvse fc Norrisl 120031) . Therefore, the ques- 
tion presents itself, is a static analytical DM potential a valid 
assumption for a study of the disruption of satellites? 



While Zhao et al. (1999) already addressed the issue 
of the evolution of a satellite galaxy in a time-varying ana- 
lytical potential, we are going to complement this study by 
comparing the disruption processes of satellite galaxies in 
analytical potentials and fully self-consistent cosmological 
simulations. 
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2 THE SIMULATIONS 

This study focuses on halo #1 from the seri es of high- 
resolu tion TV-body galaxy clusters simulations ijOill et al. I 
l2004allcl) . Halo #1 was chosen because it is the oldest of the 
simulated clusters (i.e. 8.3 Gyrs) having a reasonably quiet 
merger history. The self-consistent cosmological simulations 
were carried out using the publicly available adaptive mesh 
refinement code MLAPM jKnebe. Green fc Binnev 1200 ill 
standard ACDM cosmology (Qo = 0.3, f2 A = 0.7, Q b h 
0.04, h = 0.7, as = 0.9). For a more details and a n elaborate 
study o f these simulations we refer to reader to lGill et al. I 
j2004aO . 

From this simulation we chose one particular satellite 
galaxy orbiting within the live host halo and refer to its evo- 
lution as the "live model" . This choice was based upon two 
constraints: firstly, the satellite contains a sufficient number 
of particles and secondly, it has had multiple orbits. The 
satellite used throughout this study complies with these cri- 
teria in a way that it consists of ~ 15,000 particles (as op- 
posed to ~ 800,000 for the host) at the "initial" redshift 
z — 1.16 (8.3 Gyrs ago) and has roughly 4 orbits within the 
host's virial radius until z — when the mass of the host 
reached 2.8 x 10 14 /!," 1 M (roughly 1,550,000 particles). 

After extracting the satellite from the cosmologi- 
cal simulation we used a tree TV-body code (GCD+: 
iKawata fc Gibson I I2003T) to model its evolution for 8.3 
Gyrs in two (external) analytical potentials. The first of 
these two models we call the "evolutionary model" which 
uses an analytical reconstruction of the live halo's po- 
tential as follows: at each available snapshot of the live 
model, we fitted the host's DM density profile to the func- 
tional form of a (spherical) Nav arro, Frenk & White profile 
jNavarro. Frenk fc White1ll997l NFW): 



P(r) 
Pb 



(r/r 3 )(l + r/r s )2 



(1) 



where pb measures the cosmological background density, 8 S 
controls the amplitude and r s measures the radius where the 
profile turns from its logarithmic slope of dlogp/dlogr = 
— 1 to dlog p/dlog r = —3. From this series of snapshot 
NFW profiles we reconstruct the evolution of the param- 
eters 8s{z) and r s (z). Our best fitting functions are given 
below: 



23456 



(z + 0.08)°- : 



1015 , 



(2) 



167.28 - 58.14 z L 



Fig. shows the evolution of r s , S s and consequentially 
mass as a function of redshift for both the numerical sim- 
ulation alongside the analytical formula described by equa- 
tion 0. The increase of the host's mass around z = 0.45 
is related to a transient ( "backsplash" ) satellite galaxy 
passing just within the viria l radius at high velocity (cf. 
iGill. Knebe fc Gibson I l2004b). The second model assumes 
a static analytical host potential and hence is labeled "fixed 
model". The parameters 5 S = 21070 and r s = 97.8 adopted 
for the fixed model agree with the initial values for the evo- 
lutionary model. Note that in both analytical models the 
hosts are assumed to be spherical. Dynamical friction is not 
implemented, either. 




Figure 1. The redshift dependence for r s (upper panel), <5 S (mid- 
dle panel) and consequentially virial mass, M vlI (lower panel). 
The solid lines indicate the evolution in the live model, and the 
dashed lines are our best fitting formula to describe the evolution 
(see eg, I5t. 



3 THE RESULTS 

In Fig. |5] we show the real-space distribution (left panel) 
of the disrupted satellite after the 8.3 Gyrs evolution, i.e. 
at z — along with its distance to the host as a function 
of time (right panel). The spread in the number of orbits 
amongst the models can readily ascribed to the difference 
in host mass and the small "mis-modeling" of its growth as 
seen in Fig. at around redshift z = 0.3. 

This figure further highlights a number of interesting 
differences and similarities amongst the models. The most 
striking feature is that neither of the analytical models is 
capable of producing the real-space distribution of particles 
seen in the live model, and the live model shows the most 
"compact" distribution. Although all models display "shell" - 
like features, they are less eminent for the live host. Instead, 
the live model exhibits a "cross-like" feature and appears to 
be more compact in the central region, respectively. A com- 
parable feat ure, in fact, has a lso been noted in observational 



streams fcf. lHau et al. I 



iHelmi fc de Zeeuw 



2004). 



2000) outlined a method for iden- 
tifying stellar streams within observational data sets by util- 
ising conservation of energy, E, and angular momentum, L, 
(i.e. the integrals-of-motion) for spherically symmetric and 
time- independent potentials. In fact, this method proved 
to be a powerful tool for identifying streams in the Milky 
Way using the proper motions of solar neighbour stars (e.; 



Helmi et al. Il999l:lchiba fc Beers 120001 iBrook et al. 12003: 

Navarro. Helmi fc Freeman II2004T) . We expand their analy- 
sis not only limiting ourselves to static analytical models but 
extending the study to include both our live and evolution- 
ary models. 

Before exploring the evolution in the E — L plane (here- 
after also called "integral-space") of our target satellite in 
all three models, we present the evolution for eight different 
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Figure 3. The distribution of satellite particles in the E — L plane for the live models. The left panel shows the distributions at the time 
the respective satellite galaxy enters the virial radius of the host whereas the right panel presents the distributions at z = 0. Different 
colours represent particles of different satellites. 



satellites in the live model alone. The result can be viewed 
in Fig. |3]where each satellite is represented by an individual 
colour with our target satellite plotted as cyan. The satel- 
lites shown in Fig. [3] are all taken from the self-consistent 
cosmological simulation and represent a fair sample of all 
158 available satellites, i.e. their masses cover a range from 
from M ~ 5 X 10 10 /i _1 M to M ~ 4 x lO 12 ^ 1 M Q and 
they are on different orbits. The left panel of Fig. |3] shows 
the distributions at the time the satellite enters the virial 
radius of the host, whereas the right panel displays the dis- 
tributions at 2 = 0. 

Fig. 01 allows us to gauge the "drift" of satellites 
in integral-space. We note that, compared to Fig. 4 
of iHelmi fc de Zeeuw I l)200(f) , the integrals-of-motion are 
hardly conserved in the live model, neither for high or low 
mass streamers. The distributions rather show a large scat- 
ter, and have been significantly "re-shaped" over time. In 
addition, the mean values of E and L are also moved after 
the evolution. For instance, the "red" satellite drifts in time 
over to the initial position of the "cyan" satellite. Fig. ^fur- 
ther demonstrates that the drift of our target satellite (cyan 
dots) is comparable to the evolution of the other satellites, 
thus indicating that this target satellite is a "typical satel- 
lite" in that respect. 

One encouraging result implied from Fig. [3] however, 
is that even though the integrals-of-motion are chang- 
ing over time, satellites still appear coherent in the E — 
L plane. Hence, the inte gral-space analysis pioneered by 
IHelmi fc de Zeeuw I i2000fl still proves to be a useful diag- 
nostic to identify streams. We further like to stress that ob- 
servations only provide us with the snapshot of the distribu- 
tion at today's time, i.e. the right hand panel of Fig.|3]and 
hence measuring "evolution" is beyond the scope of RAVE 
and GAIA. 



Fig.0|now focuses on our target satellite alone, showing 
its particles at initial and final time in the E — L plane for 
all three models, i.e. the live, evolutionary and fixed model. 
We notice that particles tend to form "stripes" in the E — L 
plane (parallel to the L-axis) indicative of a spread in an- 
gular momentum for particles of comparable energy. These 
stripes form over time and have been directly linked to both 
apo- and peri-centre passages of the satellite. The live model 
deviates most prominently from the other models not only 
reducing total energies, but also "randomising" the angu- 
lar momentum and hence lacking the prominence of these 
stripes. We also observe in Fig.0|that the fixed model is the 
only model to show a noticeable number of particles (~ 3%) 
not bound to the host halo, i.e. E > 0. This can be ascribed 
to the lower mass of the host halo. 

Moreover, Fi g. H] confirms the findings of 
IHelmi fc de Zeeuw I <l200Cf) that within a fixed host 
potential a satellite galaxies retains its identity in the 
E — L plane and the integrals-of-motion hardly change, 
respectively. This statement is further strengthened by 
Fig. which presents the actual frequency distribution of 
E and L at the initial and final time. 

Not surprisingly, the live model shows the largest de- 
viations from the initial configuration (Fig. |5J. In the live 
model the distribution of E and L is broadened over time, 
and the peak of E systematically moves toward lower values, 
while the peak of L does not change dramatically for this 
particular satellite. This "drop" of E is also observed in the 
evolutionary model; it reflects the steady mass growth and 
related deepening of the host potential (cf. Fig.0. Since the 
drift in E is the most significant change, the evolutionary 
model displays a similar distribution to the live model in the 
E — L plane at z = 0. 

To further quantify the evolution of the satellite in 
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Figure 2. The projected particle distribution for all the models 
at z = 0. The lines indicate the orbital paths of the satellite of 
the respective model. 



integral-space, we calculate the centre and area of the par- 
ticle distribution seen in Fig. 0] as a function of time. The 
centre is defined to be the two dimensional arithmetic mean 
in E and L values, respectively. The area covered in the 
E — L plane is computed using a regular 128 3 grid covering 
that plane from L m i n to L max and E m - m to -E max in the re- 
spective model. Here, only cells containing more than three 
particles are taken into account. The evolution of both area 
and centre as a function of time since z = 1.16 is presented 
in Fig.© 

The top panel of Fig. |U (showing the change in area) 
demonstrates that all three halos loose coherence of the 
distribution in integral-space. This loss is most prominent 




Figure 4. Total energy, E, versus absolute value of the angular 
momentum, L, for all satellite particles within the target satellite 
at initial redshift z = 1.16 (black) and at z = (blue) in the 
live (left), evolutionary (middle) and fixed (right) models. The 
horizontal line indicates E = 0. 

for the live model and least for the fixed model. This phe- 
n omenon can also be s een in a comparison between Fig. 4 
of iHelmi fc de Zeeuw I (|200C|) and the right panel of Fig. EH 
Therefore, we conclude that a fixed (and even an evolution- 
ary) model significantly underestimates the scatter in the 
integral-space, and earlier predictions based on the simula- 
tion using a static (or time-varying) analytical description 
for the halo potential overestimate their efficiency of the de- 
tection ability of streams. 

The bottom panel of Fig. |5] (showing the drift of the 
distribution's centre in the E — L plane with respects to its 
initial position) is yet another proof that the fixed model 
holds the best conservation of energy and angular momen- 
tum. The "live" nature of the host mass can be held respon- 
sible for not only the increased dispersion (i.e. area) but 
also for the actual drift in integral-space. We further like to 
stress that at not time the evolutionary model matches the 
effects of the "live" model. There always appears to be a 
well pronounced discrepancy which can be ascribed to ef- 
fects not considered in the evolutionary model, such as the 
triaxiality of the host and interactions of the debris with 
other substructure. This leads to the immediate conclusion 
that for an understanding of the formation of streams in 
integrals-of-motion space it appears rather crucial to take 
into account the full complex evolutioary effects of the host. 
The first such step has been taken by Zhao et al. (1999) who 
modeled the evolution of tidal debris in a time- varying po- 
tential. However, our study further indicates that in a "live" 
scenario there are additional effects at work leading to an 
even greater dispersion. 



4 CONCLUSIONS 

In hierarchical structure formation scenarios as favoured 
by recent estimate s of the cosmological parameters 
JSpereel et al. 1 2003) DM halos continuously grow via both 
merger activity and steady accretion of material. Further- 
more, they also contain a great deal of substructure and 
are far from spherical symmetry. We have investigated the 
development of streams due to the tidal disruption of satel- 
lite galaxies in a DM halo forming in a fully self-consistent 
cosmological simulation. This not only models the time- 
dependency of the potential (reflecting the mass growth of 
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Figure 5. The distribution of energies (upper panels) at initial 
(black) and final (blue) time and angular momentum (lower pan- 
els) for the three models (live, evolutionary, fixed, from left to 
right). 



the host), it also accounts for other effects such as the tri- 
axiality of the host and interactions of the debris with other 
satellites. 

Our conclusions can be summarised as follows. (1) The 
distributions of (debris of) satellite galaxies in both real and 
integrals-of-motion space are sensitive to the evolution (and 
the particulars) of their host galaxy. This puts a caution on 
studies that investigate the shape of the halo based on satel- 
lite streams obtained via simulations with static DM halos. 

(2) Even in a ACDM "live" halo, satellites still appear to be 
coherent structures in the integrals-of-motion space. How- 
ever, the coherency is smeared significantly in contrast to 
predictions from simulations using static and time-varying 
host potentials, respectively. Thus, earlier studies of the de- 
tection ability of streams using a static DM pot entials (e.g. 
iHelmi fc de Zeeuw I 120001 iHarding etal~l l200ll) and even 
time-varying potentials (Zhao et al. 1999) could overesti- 
mate the efficiency of the "integrals-of-motion approach". 

(3) In the integrals-of-motion space, energy changes most 
significantly due to the mass growth of the host halo which 
is inevitable in hierarchical structure formation scenarios. 
However, there are additional effects at work such as triax- 
iality of the host and interactions of the stream with other 
satellites. Hence, any currently observed distribution of a 
satellite stream in the E — L plane no longer reflects its 
original distribution. 
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Figure 6. The evolution of the area (top panel) and the centre 
(bottom panel) in integral-space (see text for more details) as 
a function of time for the live (solid), evolutionary (dot-dashed) 
and fixed (dotted line) models. 
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